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f^ ' Abstract. We provide an overview of the Macaulay2 package VersalDeformations, which 

CSJ , algorithmicahy computes versal deformations of isolated singularities, as well as local (multi)graded 

Hilbert schemes. 
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1. Introduction 



Deformation theory provides mathematicians with the tools to describe (local) parame- 
ter spaces for various algebraic geometric objects, for example, for isolated singularities, or 

'^ ' for invertible sheaves on a projective variety, see ^HarlOj. However, computing such spaces 

J3 ■ in practice can be quite difficult, and by hand often intractable. The Macaulay2 package 

VersalDeformations, available online |Iltll] . aims to facilitate such calculations for two con- 
crete deformation problems: versal deformations of isolated singularities, and local (possibly 
multigraded) Hilbert schemes. 

\^ ■ The package VersalDeformations provides several functions which may be used to calculate 

tangent and obstruction spaces for the above-mentioned deformation problems. The function 

^ I normalMod may be used to calculate a basis for any degree of the normal module of some 

• ' (multi)homogeneous ideal in a polynomial ring. The scripted functor CT may be used to 

(^ ■ calculate bases of the first and second cotangent cohomology modules T\ and T| of some 

algebra A over a field k, assuming that these modules are finite dimensional vector spaces. 
In the (multi) homogeneous case, CT may also be used to calculate bases of homogeneous 
pieces of these modules. 

^ ■ The main contribution of the package is the method versalDe format ion, which uses 

the Massey product algorithm to iteratively lift solutions of a deformation equation to higher 
and higher order; we describe this in more detail in the following section. This can be used to 
find power series descriptions of versal Deformations and local Hilbert schemes. Since such 
a description may not be polynomial, the package provides an interface allowing the user 
to control at what point the lifting should terminate. The package also implements a more 
time-consuming lifting algorithm (via the option SmartLift) which seeks to minimize the 
number of higher order terms appearing in the equations for the parameter space. 

There are a number of other software packages which provide related functionality. J. Stevens 
has written scripts for the original Macaulay to calculate T^ and T^ for isolated singularities 
|Ste95] . There is a library for Singular by B. Martin which calculates the versal deformation 
of an isolated singularity as well as of modules |Mar99] . B. Hovinen has written a package 
for Macaulay2 which computes versal deformations of maximal Cohen-Macaulay modules on 
hypersurfaces |HovlOj . Finally, J. Boehm is developing a package for computations involving 
deformations of Stanley- Reisner rings |Boel9j . 
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2. Solving the Deformation Equation 

In the following, we briefly describe the Massey product algorithm as we have implemented 
it. For more details and mathematical background, see |Ste95j or |Ste03j . For simplicity, we 
restrict to the case of the versal deformation of an isolated singularity, although our approach 
for Hilbert schemes is similar. 

First we fix some notation. Let S* be a polynomial ring over some field k, and let / be an 
ideal of S defining a scheme X = Spec S/I with isolated singularity at the origin. Consider 
a free resolution of S/I: 

■■■ > S^ -^ S"" -^ S > S/I > 0. 



Let (pi G Hom(5"'"/ Imi?", S) i = 1, . . . ,n represent a basis of 

T^/j ^ Hom(5™/Imi?°,5)/ JacF°. 

We introduce deformation parameters ti, . . . ,tn and consider the map F^: S'[ t ]"* — )■ S[t] 
defined as 

n 

Let m be the ideal generated by ti, . . . , t„. It follows that there is a map R^ : S[t]'' — ^ •S'ft]™ 
with R^ = R^ mod m satisfying the first order deformation equation 

F^R^ = mod m^. 

Our goal is to lift the above equation to higher order, that is, for each i > 0, to find 
F': 5[t]™ -> S[t] with F' = F'-^ mod m^ and R' : S[tf -^ S[t]"" with R' = R'-^ 
mod m* satisfying F*i?* = mod m*^^. In general, there are obstructions to doing this, 
governed by the (i-dimensional k vector space T|,j. Thus, we instead aim to solve 

(1) {F'R'f + C'-^G'~^ = ^ modm^+\ 

Here, G"~^: k[t] -> k[tY ^"^^ C**~^: S[tY -^ S[tY are congruent modulo m* to G""^ and 
C*~^, respectively. Furthermore, we require that G* and C* vanish for z < 0, and C° is 
of the form V ■ D, where V G Hom(5''^, S'') gives representatives of a basis for Tg,j and 

D G Hom(S''^, S'^) is a diagonal matrix. The G* now give equations for the miniversal base 
space of X. 

Our implementation solves ([T]) step by step. Given a solution (F*, i?*, G*~^, C*~^) modulo 
m*"*"^, the package uses Macaulay2's built in matrix quotients to first solve for F*+^ and G*~^ 
(by working over the ring S[t]/I + Im(G'~^)*'' + m*"*"^) and then solve for i?*"*"^ and C*~^. 
For the actual computation, we avoid working over quotient rings involving high powers of 
m by representing the {F\ R^, G*~^, C*~^) as lists of matrices which keep track of the orders 
of the tj involved. 

3. Examples 

We provide two examples: a versal deformation and a multigraded Hilbert scheme. We 
begin with the classical example of the miniversal deformation of the cone over the rational 
normal curve of degree 4, see |Pin74] . 



VERSALDEFORMATIONS 

il : loadPackage "VersalDef ormations" ; 

12 : S=QQ[x_0. .x_4]; 

13 : I=mlnors (2 , matrix { { x_0 , x_l , x_2 , x_3 } , { x_l , x_2 , x_3, x_4 } } ) ; 

03 : Ideal of S 

14 : FO=gens I; 
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04 : Matrix S < S 
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o5 : Matrix S < S 

We see that the tangent space T^/^ of the miniversal deformation is four-dimensional, and 

■ s/i 



-s/i 
the obstruction space T|,j- is three-dimensional: 



16 : CT"1(F0; 
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o7 : Matrix S < S 

In this example, our algorithm gives a polynomial solution to the deformation equation: 

18 : (F,R, G, C) =versalDeformation (FO) ; 

Calculating first order deformations and obstruction space 

Calculating first order relations 

Starting lifting 

Order 2 

Order 3 

Solution is polynomial 

The miniversal base space is the union of A'^ and A^, intersecting in a point, as can be seen 
from the following equations: 

19 : sum G 

o9 = I t_2t_3-t_3"2 I 
I t_lt_3 I 

I t_3t_4 I 



o9 : Matrix {S[t , t , t , t ]) <- 
12 3 4 

A versal family is given by 

110 : transpose sum F 
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olO : Matrix {S[t , t , t , t ]) < {S[t , t , t , t ]) 

1234 1234 

We now consider our second example: the local description of the Hilbert scheme of the 
diagonal in P^ x P^ x P^ at the point corresponding to the unique Borel fixed ideal, see |CS10] 
for more details. 

ill : S=QQ [xl, x2, x3, yl, y2, y3, zl, z2, z3, Degrees=> 
{{1,0,0}, {1,0,0}, {1,0,0}, {0,1,0}, {0,1,0}, 
{0,1,0}, {0,0,1}, {0,0,1}, {0,0,1}}]; 



il2 : I=ideal {yl*z2, xl*z2, y2*zl, yl*zl, x2*zl, xl*zl, xl*y2, x2*yl, 
xl*yl, x2*y2*z2}; 
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ol2 



Ideal of S 



113 : (F, R, G, C) =versalDef ormatlon (gens I, 

normalModule {{0, 0, 0},gens I) , CT"2 ( { 0, 0, } , gens I ) ) ; 
Calculating first order relations 
Starting lifting 
Order 2 
Order 3 
Order 4 
Order 5 
Order 6 
Solution Is polynomial 

Note that since we were interested in the multigraded Hilbert scheme, the tangent space 
is just the degree (0, 0, 0) component of the normal module of J, and an obstruction space is 
given by the degree (0, 0, 0) component of Th^. In any case, this multigraded Hilbert scheme 
is locally cut out by 8 cubics: 
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sum G 



ol4 = I t_2t_3t_4-t_2t_4t_7-t_lt_3t_8+t_lt_7t_8+t_lt_3t_13-, 
t_lt_3t_4-t_2t_3t_4-t_lt_7t_8+t_2t_7t_8-t_lt_3t_13+. 
t_lt_3t_16-t_2t_7t_16-t_lt_14t_16+t_2t_14t_16- 
t_lt_3t_18-t_2t_7t_18-t_lt_14t_18+t_2t_14t_18- 
t_2t_4t_17-t_lt_8t_17+t_lt_13t_17-t_2t_13t_17- 
t_2t_4t_18-t_lt_8t_18+t_lt_13t_18-t_2t_13t_18- 
t_3t_4t_17-t_7t_8t_17-t_3t_13t_17+t_7t_13t_17- 
t_3t_4t_l 6-t_7t_8t_l 6-t_3t_l 3t_l 6+t_7 t_l 3t_l 6- 

There are in fact 7 irreducible components of the Hilbert scheme which pass through this 
point: 

115 : # prlmaryDecomposltlon Ideal sum G 
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